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ABSTRACT 

The effect of thrust-to-weight ratios, and apocenter altitude on 
trajectory parameters has been investigated for descent to the lunar 
surface from parking orbit. A single stage with constant thrust directed 
against the velocity vector was used in all instances. Specific impulses 
of 300 and 420 sec were used for comparison. 

The purpose of this report is to present trajectory parameters 
for the powered lunar landing maneuver from various lunar orbit 
altitudes. The ^results of the study show the variation of trajectory 
parameters for earth thrust-to-weight ratios from 0.2 to 0.5. These 
parameters include: true anomaly at initiation of powered descent; 

change in altitude, flight path angle and range angle during thrusting; 
burning time, characteristic velocity and velocity losses. 
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SUMMARY 

The effect of thrust-to -weight ratios and apocenter altitude on 
trajectory parameters has been investigated for descent to the lunar 
surface from parking orbit. A single stage with constant thrust 
directed against the velocity vector was used in all instances. Specific 
impulses of 300 and 420 sec were used for comparison. 

The results of the study show the variation of trajectory 
parameters for earth thrust-to-weight ratios from 0.2 to 0.5. These 
parameters include: true anomaly at initiation of powered descent; 
change in altitude, flight path angle, and range angle during thrusting; 
burning time, characteristic velocity, and velocity losses. 

SECTION I. INTRODUCTION 

Orbital operations are of great interest in lunar exploration 
because of requirements imposed by the landing site, energy and 
tracking. Preliminary analysis of lunar missions requires a rapid 
method of sufficient accuracy for determining trajectory parameters. 

The purpose of this report is to present trajectory parameters 
for the powered lunar landing maneuver from various lunar orbit 
altitudes. The parameters presented are: velocity, altitude change, 
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surface range, change in flight path angle, change in central angle, 
true anomaly, burning time, propellant mass fraction and velocity 
losses. 


The approach used was to apply an impulse in the parking orbit 
to attain a transfer ellipse with a 10-kilometer pericenter altitude. 

At the proper point on the transfer ellipse, powered braking was initiated 
to arrive at a 3 00 -meter altitude with the zero velocity. The equations 
of motion were integrated on a RECOMP II computer, using Runge- 
Kutta numerical integration. 

SECTION II. ANALYSIS 


One method of landing on the lunar surface is to direct the 
thrust against the velocity vector (a = 180°) and follow a gravity turn 
path. The sketch below illustrates the conditions at initiation of burning. 



First,, an impulse (as calculated from FIG 1) was applied against the 
velocity vector to attain the transfer ellipse. The pericenter altitude 
of the transfer ellipse was taken as 10 km to insure against collision 
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with the surface if powered descent was not initiated. This pericenter 
altitude was selected because of surface irregularities, guidance errors, 
retro-thrust errors, etc. From a performance view, it is preferable 
# to expend the energy at as low an altitude as possible. 

The initial true anomaly (measured from the undisturbed 
Keplerian pericenter radius) varies with thrust-to-weight ratio and 
parking orbit altitude. It is seen that as parking orbit altitude decreases 
the initial true anomaly will increase (to a maximum of 180°). 



The sketch above shows the initial and final trajectory param- 
eters, as well as change in flight path angle, altitude, and central 
range angle. The line of sight and angular range can be determined 
from FIG 2, using initial conditions. The equations of motion were 
numerically integrated from initial conditions to 300-meters above the 
surface. The initial conditions were found by iteration to give the 
desired end conditions. 

The following point mass, planar equations of motion were 
used with a = 180° in all instances. 
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V = 


_ g r 2 

F cos a €. C 


cos 
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v,/U + 


m 


g. r 


« C 


V 2 

r 


in IT 


sin 


r = V cos ^ 


• V sin 4 / 
+ = I 


(1) 

(2) 

( 3 ) 

( 4 ) 


where 


and 


m = m + mdt 
o J 


m = - 


ex 


( 5 ) 


( 6 ) 


The velocity and flight path angle may be obtained by integrating 
the equation of motion 


V = j V dt 
) 

•dt 


( 7 ) 

( 8 ) 


The range and pericenter altitude can then be calculated by 
the relations 


x = f -2* V sin if* dt 

j r 

h = h + f r dt 
o J 

and the central angle 


( 9 ) 

( 10 ) 


(ID 
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and surface range 
S = r* A4 j 

The initial weight of the vehicle is 


w = w + w 8 

o c 


( 12 ) 


(13) 


The undisturbed Keplerian true anomaly at initiation of 
powered descent is 

r \2 i 


V = COS ” 1 
O 


(V a r a ) 2 / - r c 


( 2 -"V°) sin2 A> 1 


(14) 


The velocity expended by the vehicle is the characteristic 
velocity, or 


AV = V In 
ex 




(15) 


and may be determined from FIG 3. Consequently, the velocity losses 
are the difference between the characteristic velocity and the change 
in comparative velocity, or 


AV, = AV - AV v 
loss 

where the comparative velocity is 
V" 




( 16 ) 


(17) 


and the change in comparative velocity during descent from r = r^ to 

AV ‘ = -\l V o +2|i ‘(t ' U8) 


r = r^ is 


By substituting in original equation for velocity losses 


AV, = V In 
loss ex 


(t^t )-[V v ° +2 <vi)' 


(19) 
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SECTION III. ASSUMPTIONS 

A summary of the basic assumptions used in this analysis 
follows: 


1. Deceleration of a single stage from a transfer ellipse to a 
300 meter (hover) altitude and zero velocity using constant thrust 
directed against the velocity vector. 

Z. Initial parking orbit altitudes 
h = 50 km 

h = 100 km 
h = ZOO km 
h = 300 km 

3. For comparison, specific impulses of 300 sec and 4Z0 sec 
were used. 

4. The initial earth thrust-to-weight ratio for a chemical stage 
was varied parametrically from 0.Z to 0.5. 

5. Mean spherical moon 

fj^= 4906 km 3 / sec 2 . 
r c = 1738. 3 km 

SECTION IV. RESULTS AND CONCLUSIONS 

The results are shown in FIG 4 through 13. For convenience, 
most variables are shown as a function of initial vehicle thrust-to- 
weight ratio. 

The characteristic impulse velocity, shown in FIG 4, is the 
velocity necessary to attain a transfer ellipse from a circular parking 
orbit altitude. The apocenter altitude of the transfer ellipse is 
assumed to be circular parking orbit altitude and the pericenter altitude 
is 10 km in all instances. Since losses were negligible, an impulsive 
velocity was assumed. 
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The characteristic velocity for descent to 300-meters above the 
lunar surface is shown in FIG 5a and 5b for specific impulses of 300 sec 
and 420 sec, respectively. The corresponding propellant mass fraction 
is shown in FIG 6a and 6b. 

FIG 7a and 7b show the altitude change. Ah, during powered 
descent for specific impulses'of 300 sec and 420 sec, respectively. A 
limit is reached, for a gravity turn (a = 180°), when the altitude change 
become^ equal to the apocenter altitude. This limit determines the 
minimum initial thrust-to-weight ratio for continuous burn from parking 
orbit to 300 m above the surface. This limit is reflected in the change 
in flight path angle, Atf , FIG 8a and 8b when Atf = 90°. 

The change in central angle, (angular range), during powered 
descent is given in FIG 9a and 9b. The product of A^ (in radians) and 
the mean lunar radius is the surface range and is shown in FIG 10a and 
10b. The limit shown for Ai|j and S is reached when powered descent is 
initiated at apocenter. 

The true anomaly, v Q > shown in FIG 11a and lib, is measured 
from the undisturbed Keplerian pericenter radius to the radius at 
initiation of powered descent. True anomaly equals 180° for the limit 
of minimum thrust-to-weight ratio. Burning times for the powered 
descent are given in FIG 12a, 12b, 12c, and 12d. The true gravity 
losses are shown in FIG 13a and 13b. Burning times and gravity losses 
are a maximum in the limiting case. 
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SECTION V. GRAPHIC PRESENTATION 
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Height of Spacecraft Above Lunar Surfoce (dashed arcs), km. 
Slant Range From Observer to Spacecraft (solid arcs) , km 
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FIGURE 2. CHART FOR DETERMINING ELEVATION AND SLANT RANGE 
OF LUNAR SPACECRAFT (0 to 1100 km) 



Characteristic Velocity, AV (m/sec) 
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FIGURE 3. CHARACTERISTIC VELOCITY AS A FUNCTION OF PROPELLANT 
MASS FRACTION AND SPECIFIC IMPULSE 


Propellant Moss Fraction, W s /W< 
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FIGURE 4. CHARACTERISTIC IMPULSE VELOCITY TO ATTAIN TRANSFER ELLIPSE WITH r = 1748.3 km 



Characteristic Velocity, AV (m/sec) 
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FIGURE 5a. CHARACTERISTIC VELOCITY FOR LUNAR DESCENT - I = 300 sec. 
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Earth Thrust -to-Weight Ratio, F/W 0 

FIGURE 5b. CHARACTERISTIC VELOCITY FOR LUNAR DESCENT - I = 420 sec. 
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FIGURE 6a. PROPELLANT MASS FRACTION FOR LUNAR DESCENT - I = 300 sec 
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FIGURE 7a. CHANGE IN ALTITUDE FOR LUNAR DESCENT - I = 300 sec. 
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Earth Thrust -to-Weight Ratio, F/W 0 


FIGURE 7b. CHANGE IN ALTITUDE FOR LUNAR DESCENT - I = 420 sec. 
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FIGURE 9a. CHANGE IN CENTRAL ANGLE FOR LUNAR DESCENT - I = 300 sec 
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FIGURE 9b. CHANGE IN CENTRAL ANGLE FOR LUNAR DESCENT - I = 420 sec 
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FIGURE 10a. SURFACE RANGE FOR LUNAR DESCENT - I = 300 sec. 
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FIGURE 10b. SURFACE RANGE FOR LUNAR DESCENT - I = 420 sec. 

sp 




26 











True Anomaly, v (deg) 



0.2 


0.3 


0.4 


0.5 


Earth Thrust -to-Weight Ratio, F/W 0 

FIGURE lib. TRUE ANOMALY AT INITIATION OF LUNAR DESCENT - I = 420 sec 
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FIGURE 12a. BURNING TIME FOR LUNAR DESCENT - h = 50 km 
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FIGURE 13a. VELOCITY LOSSES FOR LUNAR DESCENT - I = 300 sec. 
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